CMB Anisotropy Induced by Cosmic Strings 
on Angular Scales > 15' 
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We have computed an estimate of the angular power spectrum of the Cosmic Microwave Background 
(CMB) induced by cosmic strings on angular scales > 15', using a numerical simulation of a cosmic 
string network; and decomposed this pattern into scalar, vector, and tensor parts. We find no 
evidence for strong acoustic oscillations in the scalar anisotropy but rather a broad peak. The 
anisotropies from vector modes dominate except on very small angular scales while the tensor 
anisotropies are sub-dominant on all angular scales. The anisotropies generated after recombination 
are even more important than in adiabatic models. We expect that these qualitative features are 
robust to the varying of cosmological parameters, a study which has not yet been done. 
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The spectrum of CMB anisotropy on the angular scales 
somewhat smaller than that subtended by the horizon 
at last scattering provides a powerful probe of the na- 
ture of inhomogeneities and matter in our universe [Q. 
On these scales hydrodynamical effects can leave charac- 
teristic signatures of adiabatic or isocurvature perturba- 
tions |2|, active or passive perturbations |||, an open or 
closed universe, or even a high or low Hubble constant. 
Cosmic strings are topological defects which may have 
formed in the very early Universe and may be responsi- 
ble for the formation of large scale structure observed in 
the Universe today [Q . While cosmic string induced per- 
turbations are clearly both isocurvature and active the 
interplay between these two properties allows a range of 
possible behavior for the degree scale anisotropy || . One 
method for determining the signature for cosmic strings 
is via simulation, i.e. by numerically evolving a network 
of cosmic strings in a simulated universe and comput- 
ing the pattern of anisotropy they would produce. We 
present the results of just such a calculation in this Let- 
ter. While we have not explored the full range of cosmo- 
logical parameters, or included all the effects on network 
dynamics, our results are suggestive of what is and is not 
important for small scale anisotropies from the cosmic 
string model. 



I. METHODOLOGY 

Cosmic strings only interact gravitationally with the 
rest of the matter and produce only small metric pertur- 
bations. Thus we may ignore back-reaction of the metric 
perturbations on the strings and solve for the evolution 
of the string network in an unperturbed cosmology. This 
"stiff-source" approximation allows us to do the string 



simulation and later compute the perturbations to the 
matter, ft also allows us to reuse the the same string sim- 
ulation to compute the perturbations in cosmologies with 
different matter content. Here we have reused one of the 
simulations that was was used in Allen et al. ||. While 
that simulation took the size of the simulation box to be 
twice the present day horizon, we may use the assumption 
of self-similar evolution of the string network ( "scaling" ) 
to rescale the box to a smaller size while simultaneously 
decreasing both the starting and ending time of the simu- 
lation. The equations of motion of the strings will still be 
satisfied so long as the cosmological expansion remains a 
power law, which in this case is a cx fs , corresponding to 
a flat FRW matter dominated universe. These simula- 
tions should not be used for epochs too close to matter- 
radiation equality, which means that if we want to study 
the effects close to recombination we should really only 
consider a large Hubble constant which puts recombina- 
tion long after matter-radiation equality. Here we use 
Hq = 80, but even with this Hubble constant the effects 
of the matter-radiation transition are liable to be signif- 
icant, and below we discuss what these effects are liable 
to be. Throughout we take the baryon fraction = .02 
to be consistent with the predictions of nucleosynthesis. 

The cosmic string simulation is in a cubical box whose 
comoving size we denote by L. This simulation starts 
when the comoving horizon (= conformal time) given by 
77; = 0.05L and ends at r){ = 0.5L. During this period the 
universe has grown by a factor of 100. The final time was 
chosen to the maximal value for which one is guaranteed 
that the network properties will not be affected by the 
periodic boundary conditions put on the string network. 
The statistical properties of CMB anisotropy generated 
in this box will also not be affected by the periodicity so 
long as the anisotropies are measured at conformal time 



1 



??obs < fji- I n order to use this simulation to compute 
small-scale anisotropics produced near recombination one 
must make the measurement of the anisotropies long af- 
ter the end of the string simulation. One may do this 
by rescaling the simulation to end just after recombina- 
tion, but continuing to propagate the brightness pattern 
unperturbed through the periodic box until today. This 
free- streaming procedure will lead to a brightness pattern 
which is quasi-periodic on the sky, the size of the peri- 
odic patches being the angle subtended by length L when 
placed at the distance of the surface-of-last-scattering. 
We expect that the statistical properties of the pattern 
in each of these patches will be nearly the same as for a 
similar sized patch in an infinite non-periodic universe. 

Since the metric perturbations are small we may use 
the linearized Einstein equations for the metric and the 
linearized Boltzmann equation for the anisotropies. This 
makes the equations for the temperature anisotropy lin- 
ear in the stress-energy of the strings, 8 M „. The dynam- 
ics of the strings is described by non-linear equations, 
and leads to the non-Gaussian distribution for ®uvi but 
the response of the photons to the strings is linear. The 
solution of the linear equations may be written as a ho- 
mogeneous and inhomogeneous part 

AT „ , /AT 

— (n,x obs ,?7 obs ) = I — 

J drl J dVZ^(n,x,<77,V)e^(x'n') (1.1) 

where depends on the initial condition for the cos- 

mological fluids (CDM, baryons, photons, and neutrinos) 
which is often referred to as compensation. We have com- 
pensated the strings with adiabatic perturbations in the 
other matter so as to make the total density initially uni- 
form, just as in Ref. Q. Clearly this is not exactly the 
correct prescription except on supcrhorizon scales. Nev- 
ertheless we find that (^-) l is small in comparison with 
the 2nd term (this is consistent with the results of Ref. 
j?J). We think this is likely to remain true with different 
prescriptions for compensation. 



The way we have used Eq. (1.1) is somewhat different 
than has typically been done in the past. Rather than 
concentrating on a few x Q b s we will effectively compute 
the entire brightness pattern throughout the simulation 
box. This is a function of the 5-d phase space (n, x t, s ) 
and is quite a large object. One can reduce this size 
by Fourier decomposing the x-dependence and then, for 
each k-mode, expanding the n-dependence in spherical 
harmonics, Y (l m) , using a spherical polar coordinate sys- 
tem with the pole in the direction of k. The beauty of 
this decomposition is that the m — terms correspond 
to scalar modes, the m = ±1 terms to vector modes, and 
the m = ±2 to tensor modes ||. Modes with |m| > 2 
are zero since they are not coupled to the gravitational 
field (in linear theory) and may be ignored. From the 
mode coefficients, A (! _ m) (k), one may construct the vol- 
ume average of the C\ 's measured by all the observers in 



the box: 
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There are no cross-terms between different terms because 
all of the modes are orthogonal. Applying this decom- 
position to the right-hand-side of Eq. (1.1) induces a 



scalar/vector/tensor decomposition of D^ u and MI/ . We 
may perform the decomposition of M „ numerically. The 
corresponding components of are solutions of the 
linearized Einstcin-Boltzmann equation which we have 
computed numerically using standard techniques H . 



II. TESTS 

Our calculation has three parts: the string simulation, 
the numerical calculation of the Green functions, and 
the computation of the brightness perturbation, which 
merges the two preceding parts. A number of tests have 
been made on the simulation P,[l0[|7|] which we do not 
discuss here. 

The Green functions have been tested in two ways. 
First, the code which computes the scalar and tensor 
Green functions ^] has been used to determine the an- 
gular power spectrum for a model with adiabatic pertur- 
bations only. These calculations agree to within a percent 
with standard results. Although the vector Green func- 
tions cannot be tested in this way, there is an analytic ap- 
proximation in terms of spherical Bessel functions, which 
holds after recombination. For this case, there is excellent 
agreement between the analytic and numerical results. 
Similarly, the scalar Green functions have been tested on 
large scales by comparing with analytic approximations; 
again we find excellent agreement. 

To test the full pipeline of merging the Green functions 
with the stress energies, we calculated the anisotropies in 
a matter dominated universe with one domain wall [fll| . 
The C/'s in this case can be calculated analytically. The 
numerical results are in very good agreement with the 
exact results. 

The greatest limitation of our ability to accurately 
compute the small angle CMB anisotropy is the lack of 
dynamic range of the cosmic string simulation. Although 
the simulation runs over a range of 10 in conformal time, 
this is insufficient to accurately follow the evolution of 
modes from well outside to inside the horizon, where they 
begin to oscillate. Our results indicate that modes well 
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inside the horizon make a significant contribution to the 
anisotropy. 

Due to the limited spatial resolution of the cosmic 
string simulation, a limited number of Fourier modes 
of the source stress energy are available. To determine 
whether this would prevent us from seeing "Doppler 
peaks," we used the same number of modes over the same 
range of time for a model with adiabatic initial condi- 
tions, for which the exact results are known. Qualitative 
agreement was obtained, as the peak structure and loca- 
tion were clearly evident. Hence, we do not expect that 
the limited spatial resolution should prevent us from re- 
solving features in the angular power spectrum. Finite 
grid effects lead to an artificial drop of only ~ 10% in the 
amplitude of the angular power spectrum for the largest 
Z-values obtained from a simulation volume. 

Finally, we have used the same simulation to calculate 
anisotropies generated at different times. However, we 
do not expect there to be strong correlations in cither 
the source stress energy or the anisotropy on time inter- 
vals longer than the run time of the string simulation. 
Hence, we do not expect the re-cycling of the numerical 
simulation to have a strong effect on our results. 

III. RESULTS 

We obtain the final angular power spectrum, C;, by 
combining the results from different rescalings of the re- 
sults of the cosmic string numerical simulation. As dis- 
cussed earlier, we exploit the scaling properties of the 
string network to make up for the finite temporal dy- 
namic range of the simulation. For multipole moments 
in the range 2 < I < 20 (large angles) we use one box, cov- 
ering the redshift range < z < 100. For 100 < I < 800 
(small angles) we use two different boxes corresponding 
to 700 < z < 17000 and 60 < z < 700. 

Our results are shown in Figures 1 and 2. On large 
angular scales, we see that the dominant source of the 
anisotropy is due to vector perturbations. This is an 
important result, as it means the COBE normalization 
of the mass-per-unit-length fixes the amplitude of the 
vector, not scalar, spectrum. 

We may use the large angle results as a test of our 
techniques, by comparing with the results of Ref. Q, 
which used the same simulation. They found i(i + l)C; ~ 
350(G^) 2 at large angles, so our results are a factor of 
~ 2.8 smaller. However, we are in accord with the shape 
of the spectrum and the determination that the major- 
ity of the large angle anisotropy is generated at redshifts 
z < 20 (as in Fig. 3 of Ref. ||). When normalizing 
to COBE fl|, we find ^ 6 = Gfi x 10 6 = 1.7, a fac- 
tor ~ V2.8 higher than that of Ref. Q. This result is 
comparable with ^ = 1.5(±0.5) flf], 1.7(±0.7) @ and 
2 p!q |. Thus, our large angle normalization falls in the 
middle of the range spanned by other calculations. 



On small angular scales Figure 1 shows the contribu- 
tion to the scalar Ci's from early times (z > 700). There 
is a gradual rise in the spectrum from Z = 100 till a very 
broad plateau for I > 200. The increase is less than a 
factor of two so may not be significant. 

Calculating the contribution from the vector modes 
and from the late time contributions to the scalar modes 
presents a problem. In each case, the box size must be 
chosen very large in order to get to very late times. The 
smallest scales (I > 300) therefore become unreliable. 
In Figure 2, we show the results for 100 < I < 300. 
The most striking feature here is that both the vector 
contribution and the late time scalar contribution are 
larger than the early time scalars. The net result is that 
the mild acoustic peak in the scalar spectrum becomes 
hidden. The total spectrum therefore is quite flat. (It is 
unclear if the slight drop in the total is significant.) While 
these qualitative features appear to be robust, we believe 
that in future work they will be subject to quantitative 
changes which will tend to significantly boost power on 
small- angular scales. 

IV. DISCUSSION 

We have calculated the anisotropies in the CMB and 
the perturbations to matter induced by a network of cos- 
mic strings in the matter era. This calculation uses all the 
components of the stress-energy tensor and exact Green 
functions for all modes. The greatest source of uncer- 
tainty in our results is due to the limited spatial resolu- 
tion and run time of the numerical simulation. 

On large scales, we find that the vector perturbations 
are very important for the CMB anisotropy from cos- 
mic strings. Our results for anisotropies at small scales 
appear at present to be too low to be consistent with cur- 
rent measurements. Future work will determine whether 
inclusion of a radiation-matter transition in the cosmic 
string simulation, in which the long string density and 
rms velocity relaxes from the radiation- to matter-era 
scaling values, inclusion of a "wiggly" equation of state 
for the cosmic strings, or a lower value of the Hubble 
parameter and a larger value of the baryon density can 
boost the small-scale power sufficiently. Even if a pri- 
mary Doppler peak were to become a more prominent 
feature, this work suggests that secondary Doppler peaks 
are not anticipated in cosmic string models. Finally let us 
mention that qualitatively similar results have also been 
found for global defects Eg], 
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FIG. 1. Angular Power Spectrum at large and small scales 
for the cosmic string simulation. Tensors do not contribute 
appreciably. On small scales only the contribution from early 
times (z > 700) can be calculated reliably over the full range 
of angular scales. 



FIG. 2. Angular Power Spectrum at small scales for the 
cosmic string simulation. The three different contributions 
are: the early time scalars; late time scalars; and vectors. 
These add incoherently to give the total spectrum. 
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